Conductance beyond the Landauer limit and charge pumping in quantum wires 
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Periodically driven systems, which can be described by Floquet theory, have been proposed to 
show characteristic behavior that is distinct from static Hamiltonians. Floquet theory proposes to 
describe such periodically driven systems in terms of states that are indexed by a photon number 
in addition to the usual Hilbert space of the system. We propose a way to measure directly this 
additional Floquet degree of freedom by the measurement of the DC conductance of a single channel 
quantum point contact. Specifically, we show that a single channel wire augmented with a grating 
structure when irradiated with microwave radiation can show a DC conductance above the limit of 
one conductance quantum set by the Landauer formula. Another interesting feature of the proposed 
system is that being non-adiabatic in character, it can be used to pump a strong gate-voltage 
dependent photo-current even with linearly polarized radiation. 
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I. INTRODUCTION 

Harmonically time-varying Hamiltonians offer an in- 
teresting range of phenomena such as the coherent de- 
struction of tunnelingfi photon-assisted tunneling;^ adi- 
abatic pumping^^— and non-adiabatic pumpingi^ Re- 
cently it has been realized that such periodically time- 
dependent Hamiltonians can possibly be used to real- 
ize topological phases in materials, which were origi- 
nally non-topologicalii"— Systems with periodically time- 
varying Hamiltonians (i.e. H{t + T) = H{t)) are de- 
scribed in terms of a set of quasi-stationary quasi-energy 
eigenstates e'^^*ip{t), which unlike stationary states of 
static Hamiltonians vary periodically in time ip{t + T) = 
ipit) and have energy eigenvalues E that are only de- 
fined modulo J7 = 2tt/T, where T is the period of the 
perturbation^^ These quasi-energy eigenstates, which we 
will refer to as Floquet states, can be determined from a 
time-dependent Schrodinger equation 



{H{t) - idtMt) = E^{t), 



(1) 



where the quasi-energy eigenvalue E satisfies the condi- 
tion E e [-§,§]. 

As a special case we may consider a time-independent 
Hamiltonian H{t) — H{t = 0). In this case, the fre- 
quency fl is arbitrary, and the Floquet states are just 
eigenstates of the Hamiltonian. If Es is the correspond- 
ing eigenvalue of then the quasi-energy Eg is given 
by Eg = Es — ngJl, where is the unique Floquet-index 
integer needed to ensure that Eg G [— Note that 
we have considered a time-independent Hamiltonian for 
illustrative purposes only. This case is slightly patholog- 
ical since the frequency Q, is not uniquely defined for a 
time-independent Hamiltonian, and in principle one can 
define the quasi-energy Es relative to any frequency 51 
one chooses. On the other hand, the Floquet index rij, 
defined in terms of the energy eigenvalues does not 
apply to general periodically time-dependent Hamiltoni- 
ans where the energy is not well-defined. However, 
the effectively stationary states, which are classified by 



an energy e in the time-independent case, must be clas- 
sified by quasi-energy in a limited range E € in 
the general periodically time-dependent case. Therefore 
in general, one expects a much larger density of Floquet- 
eigenstates, which can participate in processes associated 
with a specific quasi-energy eigenvalue E. The increased 
density of states at a fixed quasi-energy E can be thought 
of as resulting from the fact that a system in a putative 
state with energy e can be driven into a state with energy 
E+nVL and back by a strong time-dependent Hamiltonian. 
The expanded state space accessible to states with a fixed 
quasi-energy E, which results from the possibility of co- 
herent absorption and emissions of quanta with energy 
57, is what is responsible for the rich phenomenology of 
periodically driven systems. 

States characterized by fixed energy or quasi-energy 
play a crucial role in scattering and transport processes. 
Since the density of states at a fixed quasi-energy E in 
a periodically driven system is dramatically increased as 
compared to static systems, one expects the transport 
property of driven systems to be significantly different 
from static systems. Such dramatic differences have al- 
ready been seen in micro-wave resistance oscillation ex- 
periments in two-dimensional electron gases in magnetic 
fields. ^^'^^ More specifically, if one considers the disper- 
sion of a single-channel quantum wire in the Floquet rep- 
resentation, which is written as En^k — £k — nO., one 
immediately observes that the number of Floquet chan- 
nels of the wire is significantly larger in the 'quasi-energy' 
Floquet representation as compared to the static repre- 
sentation with fixed energy eu- 

A natural question to ask is whether the increased 
effective number of Floquet channels in a periodically 
driven one-dimensional quantum system modifies the 
transport in a direct way. In fact, the Landauer formula 
for coherent electronic transport through a static quan- 
tum systemji^^— relates the conductivity to the num- 
ber of channels in a transport system. In particular, 
the conductance of an N channel quantum wire of non- 
interacting electrons is bounded above by NGq, where 
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Go = 2e^//i is the quantum of conductance. The bound 
on the conductance is reached for quantum wires which 
are non-interacting and are connected to reservoirs by 
perfectly transmitting contacts. The increased number 
of Floquet channels for driven systems suggests that the 
Floquet channels could directly manifest themselves by 
contributing to the conductance leading to a conduc- 
tance excess of Go in driven one-dimensional systems, 
at least in principle. In this paper, we show that this is 
indeed the case, and a single channel periodically driven 
quantum wire with appropriately designed contacts can 
have a DC conductance that is in excess of the Lan- 
dauer bound Go. Moreover, for geometries which are 
not symmetric between the left and right leads, the ap- 
plication of a time-periodic drive will be found to result 
in a DC pumping current even at zero voltage bias, sim- 
ilar to previously studied systemsj^i^ Since our system 
will need to be driven at a finite frequency, we will find 
that even a quantum wire driven by linearly polarized 
microwave radiation displays a gate-voltage dependent 
pumping current. This is in contrast to the adiabatic 
pumping in the limit of small driving frequency where 
it has been shown that a two-parameter drive such as 
one resulting from circularly polarized radiation is neces- 
sary to drive a significant pumping current. — A similar 
enhancement of the conductance in strongly-interacting 
driven quantum wires has been predicted using the 
Bosonization approach. Moreover, it has been suggested 
that such an enhancement might apply to the weakly 
interaction case.— Also, single-parameter non-adiabatic 
pumping has also been demonstrated in experiments, al- 
though it is attributed to a mechanism different from the 
one discussed in the present paper.^® This paper provides 
a unified picture that is based on the Floquet-state in- 
terpretation for enhanced conductance for weakly inter- 
acting systems and correspondingly strongly gate voltage 
dependent pumping. 



Outline 

We start by reviewing in Sec. II the scattering matrix 
formalism in Floquet space. Following this in Sec. Ill, 
we discuss the analog of the Landauer formula and the 
resulting conductances for a two-terminal geometry in 
a driven system4iii^ In this section, we also review how 
unitarity of the scattering matrices leads to the bound 
on the conductance obtained by Landauer and show how 
in principle this bound may be violated in systems with 
time-dependent Hamiltonians. In Sees. IV and V we 
discuss the details of a specific wire geometry that is cal- 
culated to support, in principle, a conductance in excess 
of the Landauer limit. Furthermore, the same structure 
is found to carry a pumping current. We estimate pa- 
rameters for a realization of this structure in a GaAs 
two-dimensional electron gas (2DEG) and show that the 
excess conductance and pumping current can be mea- 
sured for experimentally reasonable parameters. In Sec. 



VI, we compare the pumping current obtained in our 
structure to the adiabatic pumping proposals using cir- 
cularly, or elliptically polarized radiation. 

II. FLOQUET SCATTERING THEORY 

Let us now discuss the problem of scattering of elec- 
trons by a harmonically time-varying potential V, with 
frequency fl, that is localized in a region in space C 
(shown in Fig. [T]) separating a left-lead L and a right 
lead R. Each of the leads are assumed to support multi- 
ple propagating channels indexed by p. The periodically 
time-dependent Hamiltonian for the multi-channel sys- 
tem is written as 

Hp'At) + V - ^F)Sp',p + [yp',p(x)e'"* + h.c], 

(2) 

where x is the coordinate along the wire, p\p are channel 
indices, m is the effective mass of the electrons, Ap is the 
confinement energy induced mini-gap between the vari- 
ous channels, ep is the Fermi-energy of the electrons in 
the wire and Vp'^p{x) are the position-dependent matrix 
elements of the potential applied in the central region. 

The solutions of the corresponding Schrodinger equa- 
tion, Eq. [1] can be written in terms of the time- 
independent wave-functions (pn.pix) 

V'p(x,t)=^(/.„,p(x)e^""*, (3) 

n 

where the index n is referred to as the Floquet index. 
The wave- functions (j)n,p{x), are are solutions of the time- 
independent Floquet equation,^^ which is written as 

T^^KZ^^^ - nn)cj,n.,p' - Ecj^^^p, (4) 

n,p' 

where H{t) = 7j("-)g™sit^ p^j. ^j^g specific case of 
the quasi one-dimensional wire described by Eq. [2l the 
components T?*^"-* are written as 

^(0) _ ( \ X , 

ij(i) = i7(-i)t = Vp,,p{x). (5) 

Since the time-dependent potential V vanishes outside 
the central region G, the solutions in the leads to Eq. |4] 
can be expanded in terms of plane-wave energy eigen- 
states 

^nA^) = e'''"''"0n.p(O) (6) 

at a pre-determined quasi-energy E, where the wave- 
vectors kn^p are determined by the requirement that the 
Floquet energy £p{knA — £f — nfl = E, with ep{k) = 
^ — Ap being the dispersion of channel p in the wire. 
Here E G [—$7/2, 17/2] is the Floquet energy of the scat- 
tering states and is the Fermi-energy of the incoming 
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FIG. 1: Schematic of conductance set-up from a 2 channel 
left lead L to a single-channel right lead R through a 2 channel 
intermediate region C. The transmission matrix To is static 
such that one of the 2 channels in C is perfectly transmitted 
into R while the other channel in C is reflected back into C. 
The floquet matrix TF(r2) contains a harmonic frequency but 
is refiectionless and can only transfer electrons between one 
channel and the other. 

electrons in the lead in the absence of a bias voltage. 
Note that the absence of a time-dependent potential in 
the leads is crucial in allowing us to consider states with 
a definite Floquet index n. Thus in the presence of a pe- 
riodically time-dependent scattering, in addition to the 
channel index p, one must label incoming states and out- 
going states in the lead with a combination {n,p) of the 
Floquet index n and the usual channel index p. The sign 
of the wave-vectors kn^p in each of the leads L, R deter- 
mines whether the mode is moving into or out of the 
junction region C. Since e^? is the Fermi-energy of the 
incoming electrons in the leads, the Floquet index n for 
the incoming modes must satisfy n = 0. We will denote 
the Floquet index of the out- going state by n'. 

Thus, the solutions to Eq. 2] can be represented con- 
veniently in terms of a scattering matrix (S-matrix) 
Sn',p';n=o,p of statcs from the incoming modes with Flo- 
quet index n = in channel p to out-going modes with n' 
in channel p' at a fixed floquet energy E. The S-matrix 
is defined by the equation 

«n>''/'iv^ = ^n'p',n=0.p{^^ E)v]l2o ,p4>n=0 ,p (7) 
P 

where w„=o,p are the group velocities in the respective 
channels. The scattering matrix S is unitary as a conse- 
quence of particle-number conservation. 



III. UNITARITY OF THE 5-MATRIX AND THE 
LANDAUER FORMULA 

The current flowing I from the left-lead L to the right 
lead R in the geometry in Fig. [1] is a function of the 
Fermi-energy e^? of the incoming modes of the leads and 
the bias voltage Vb- In the presence of a bias voltage, 
Vf), the Fermi-energy of the incoming electrons in the 
left lead, eF,L, differs from the Fermi-energy of the right 
lead €f,r. = £f,l — Vb- The conductance a{eF), at an 
equilibrium Fermi-energy 6^ is defined as the response of 
the current /(eF,I4) to an infinitesimal bias voltage Vh 
i.e. a{eF) = dI{eF,Vb)/dVb\vh=o- Therefore the current 



I{^F, Vb) is a function of the mean Fermi- level ep and the 
bias voltage Vb- 

However, to derive the Landauer formula it is conve- 
nient to consider the current / as being a function of the 
Fermi-level of the incoming electrons in the left and right 
lead which are written as 

CF.L = eF + Vb/2 (8) 
CF.fl ^ep- Vb/2. (9) 

Non-interacting electrons coming from the left lead L 
remain decoupled from electrons coming in from the 
right lead R. Therefore, the current I{cf,Lt^f,r) can be 
separated as I{eF,L,£F,R) = iLi^F.L) + iRi^F.R) where 
Il{^f,l) {Ir{^f.r)) is the current carried by states from 
the left (right) lead L (R) with Fermi-level eF,L( ^f,r)- 

In the model for non-interacting electrons considered 
explicitly here, changing the Fermi-level in either lead is 
equivalent to applying a voltage to that lead. In this case, 
the mean Fermi-level in the wire can be easily changed 
by applying a voltage symmetrically (i.e. 5eF,L = ^cf^r) 
to the two leads. For interacting electrons, however, it 
is difficult to change the net charge density in the wire, 
so a change in the average voltage of the leads produces 
primarily a change in the electrostatic potential of the 
wire and only a slight change in the mean Fermi-energy 
ep. However, if the sample has a side gate in the vicinity 
of the wire, which will be needed in our proposal to pro- 
duce the confinement needed to create the wire in any 
case, it may be possible to produce a significant change 
in the electron density in the wire, and thus vary ep, by 
applying a large voltage difference between the gate and 
the leads to the wire. Although we consider explicitly 
only the non-interacting case, we shall write our final re- 
sults in terms of cf and Vb, rather than eF.fi and cf.l in 
the expectation that they will also be applicable at least 
qualitatively, for an interacting system. 

While the S-matrix representation is convenient be- 
cause it leads to a simple condition for probability con- 
servation, for transport between the left and right leads, 
it is more convenient to think in terms of the transmission 
or 7~-matrix, which is the part of the S-matrix defined as 

/ T^(L^L) ^(L^R) \ 

Q I n'p';np n'p';np \ f^n\ 

On'p';np— I IR^L) ^(R^R) I ' \^'^) 

\ 'n'p'\np n' p' \np J 

where TZ is the reflection-matrix. The superscripts L and 
R denote whether we have modes going from left to right 
or vice-versa. Applying the standard argument used in 
the derivation of the Landauer formula to the incoming 
states in the left-lead with Floquet energies in the range 
El — ^F,L — ef S [0,14/2], the part of the transmis- 
sion matrix Tl^r. can be used to calculate the current 
Il{^f,l) and therefore also the derivative of the dc cur- 
rent / with respect to the left lead chemical potential 
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epL, as 



dl 



dh 
de 



F,L 



(11) 



p' ,n' .p 



where the channels p summed over for the occupied in- 
coming states in L and n' and p' are, respectively, the Flo- 
quet and channel indices of the outgoing modes in R. The 
response of the current / with respect to the right chem- 
ical potential eF,R, which is defined as (7r{E) = —3^7^, 
can be calculated using an expression similar to Eq. 1111 
Note the negative sign in the definition of aR{E) accounts 
for the direction of the current / being from L to R. Since 
the current / in static systems responds only to the dif- 
ference in the chemical potential Vb — eF,L — £f,r, it is 
convenient to define the conductance across the leads as 
the response of the current to the bias voltage as 



dl 
dVb 



O-RieF) + CTLieF) 



(12) 



where is held fixed when taking the derivative. 

Static systems, in the absence of a bias voltage Vb — 0, 
are in equilibrium so that the total charge in the right 
lead R must be time- independent. Therefore the current 
I{ep,Vb = 0) vanishes in the absence of a bias voltage. 
However the driven point contact, in addition to carrying 
a current in response to a bias voltage Vb, may also car- 
ries a pumping current. The derivative of the pumping 
current is calculated as 



dip 
dep 



dh 
de 



dip a-L{ep) - apiep) 



F,L 



de 



F.R. 



(13) 

where Vb is held at = when taking the derivative. 
The absence of a pumping current Ip = in the static 
case leads to the constraint 



o-L(e_F) = a-piep) for Vt = Q. 

The pumping current Ip{ep) — I{ep,Vb 
lated using the relation 



Ip{ep) = / de<7p{e), 



(14) 

0) is calcu- 
(15) 



where e_F is the Fermi-energy at which the current is mea- 
sured. 

The above discussion has been restricted to zero tem- 
perature. Finite temperature modifications of the above 
results are obtained by averaging the above conductances 
with the function '^fi'^+^^'^^^'^^T) ^^^^^ j^^^. ^^^^^ ^-^^ 

Fermi function with temperature T and Fermi-energy ep. 
For the rest of this paper we will assume that the temper- 
ature T is significantly smaller than the driving frequency 
so that it does not qualitatively affect our results. 

The expression Eq. [TT] is a generalization of the usual 
Landauer formula, which is obtained from Eq. 1111 by 



considering the static limit where the time-dependent 
potential Vp^q{x) — vanishes, so that the Floquet- 
index n is a conserved quantity and the Floquet index 
n' — Q for the out-going states. In this limit, all incom- 
ing modes from the left lead L with energy in the range 
^F,R < {E + ep) < ep^L are occupied while all the outgo- 
ing modes (in both L and R) are occupied according to 
the scattering matrix T. Since the matrix 7~'^~*^' is only 
part of the unitary S-matrix, the sum over all outgoing 
modes (p') of 



E |r;^"'''(f^-0,£;)p< E \Sp'p{n = Q,E)\^^Nn 

(16) 

where Np is the number of channels p' in the right lead 
R. It follows that the conductance from the left lead is 
less than (Jl{^f) < Np. In particular, in the case where 
the right lead R has a single channel i.e. Np — 1, the 
conductance is bounded by (JL{ep) < 1. Using Eq. [HI it 
follows that (TR{ep) < 1 and hence cr{ep) < 1. 

In the case of a time-dependent perturbation, multiple 
values of the out-going Floquet-index n' are allowed and 
at least, in principle, the conductance can exceed 1. Such 
a conductance exceeding unity provides a direct test of 
the accessibility of Floquet modes indexed by the pho- 
ton number n. Below we discuss one scenario where the 
conductance can exceed 1. Furthermore, the geometry 
we will discuss will also have a non- vanishing pumping 
current, even for single parameter pumping. 



IV. SCHEMATIC SYSTEM 

Let us now discuss a specific case of the geometry 
shown in Fig. 1 where the conductance can exceed 
the conductance quantum. In particular, we assume 
that the portions L and C of the wire have two trans- 
mitting channels each, while the region R has a sin- 
gle transmitting channel, so that the static conductiv- 
ity is bounded above by unity. For simplicity, the scat- 
tering potential in the central region C is split into a 
time-dependent but reflection-less scattering part with 
a transmission matrix Tf{^,E) and a static, but only 
partially transmitting, part with a transmission matrix 
To- The reflection- less potential with transmission ma- 
trix Tp{^,E) can be approximately realized by a grat- 
ing potential, which only has wave-vector components 
on the scale of ±(fcF,i — fcF,2)- The transmission matrix 
To is realized in the usual way, constricting the quasi- 
one-dimensional channel adiabatically from a width that 
supports two channels at the Fermi level to a width that 
supports only one. The time-dependent part of the scat- 
tering Tp{i^,E), being almost purely transmitting, only 
transfers electrons between channels p = 1, 2 in the lead 
L and central region C. The static scattering To com- 
pletely reflects in channel p = 2 and transmits channel 
p — 1 from C to R, so that the using the Floquet Lan- 
dauer formula Eq. [11] the left conductance from L to R 
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\TF,n'l;n=Opi^,E — 0)|^ 



(17) 



0=1, 2;n' 



Using the unitarity of the Tf matrix in Eq. [171 the 

term ^j^, \Tf n'i-oi{^, E = 0)p can be replaced by 

En' \TF,n'lMn,'E = 0)|2 = 1 - E„' \TF,n'2M^> ^ = 

0)p, so that aL is re- written as 
<JL{eF) = 1 

+ \TF,n'lM^,E = 0)P - |rF,n'2;0l(f^, E = 0)|2. 

n' 

(18) 

The conductance from the right lead is entirely de- 
termined by the fact that To transmits channel 1 in 
the left lead perfectly into the right lead and reflects 
channel 2. This corresponds to a (TlIef) ~ 1 for just 
the transmission matrix Tq. Since To is the result of a 
time-independent scattering, using Eq. 1141 one obtains 
"■/{(ef) = (^l{^f) = 1, so that the only incident channel 
i.e. channel p = 1 must be completely transmitted from 
right to left. The potential for the time-dependent scat- 
tering i.e. Tf is reflectionless and therefore must also 
completely transmit these electrons in channel p — \. 
The combination of these arguments lead to the result 
that, even in the presence of an AC perturbation. 



cri?,(eF) 



1. 



(19) 



Such configurations are characterized by a derivative of 
the pumping current up satisfying 



o-p{eF) = cr(eF) - 1, 



(20) 



so that the existence of the pumping current is related to 
having a conductance a above the Landauer limit a — \. 

To maximize the conductance (j{eF) in Eq. I12[ 
one must maximize aLi^p)- Since the sum 



En'lTF.n'vM^^E 



0)1 < 1, the conductance 



o'hi^F) in Eq. \TE\ is bounded by unitarity of Tf 
by 2. Therefore, if one can construct a sys- 
tem where En' ITf,™'1;02(^^, = 0)\^ ~ 1 and 
\TF,n' 2;0ii^ : E = 0)| « then CTi(ei?) can reach 
it's maximum value of 2. If aLiep) reaches its maximum 
value of 2, the conductance cr(ei?) would reach a value 
of 1.5, which is in excess of the Landauer limit of cr = 1 
for a single channel wire. In addition, there would be a 
pumping current from the fact the existence of a non- 
vanishing differential pumping current aplep) ^ 0.5. 
Below, we discuss a realization of Tf which would 
produce such a property. 



V. CALCULATION OF Tf{^,E = 0) 

We now discuss a specific model for the time-dependent 
scattering Tf{^,E), which is shown in Fig. [H that can 



lead to the ideal amplitude scattering described in the 
previous paragraph. The transmission matrix Tplfl^E) 
can be obtained by solving Eq. 2] written as 



92 

+ A 



P 



- {nfl + E)(j)n,pi 

(21) 



where we have used the Hamiltonian components in 
Eq. [5] As mentioned previously, we will choose a poten- 
tial Vpi ^p{x) for the transmission amplitude Tf(51, E), so 
that the Fourier components at the back-scattering wave- 
vectors 2fci?_p and kp.i + fcF,2 can be neglected. There- 
fore, for frequencies and temperatures much lower than 
the Fermi energy for each of the channels, we can re- 
strict our attention to states near the right moving Fermi 
points in the two channels with Fermi-wave- vectors fcF,i,2 
in the lead L. For such states of the Hamiltonian in 
Eq. [5l it is convenient to transform the wave- functions as 
4'n,p{x) = e^''"--''^(j>n,p{x), and ignore second derivatives 
of the slowly varying (compared to the absolute Fermi- 
wave-vector kF,i,2) wave-function 4>n,p{x) so that Eg. 1211 
simplifies to 



dx4'n,p' — ~''"F,p'' 

^[V,p(x)e*'="-i-^0„_i,p + i;*,,p(a;)e''="+i-^0„+i,p] 
p 

(22) 

where — {kF,p + n-^^). Propagating these equations 

with initial values at a; = so that (pn.pix = 0) is non- 
zero only for the incoming modes with n = 0, one can 
calculate the coefficients TF.n'p'-Op{^) through the equa- 
tion 

y^^TF,n'p'-Op{^)^VF,p<Po,p{x = 0) = ^VF,p'<l>n',p'{x = L) . 
P 

(23) 

Calculating the scattering matrix Tf for a weak scat- 
tering potential V shown in Fig. [2] using the first order 
Born approximation yields insight into the mechanism 
for obtaining a large Floquet conductance. Given that 
the incident states are at photon number n = 0, the final 
states in region C, following the time-dependent scatter- 
ing, are written as 

4>i.p. [L) = ^ivp\, / dxe-'"=i-'^ V.p(2^)e''°""^o,p(0) 
To p 

4>^iy{L) = -tvp], /'^dxe-^'^--- Vl/;,,p(a;)e''="--0o,p(O), 
To p 

(24) 

for p',p — 1,2. Within the first order Born approxima- 
tion, only the scattering processes shown by arrows in 
Fig. [21a), which absorb or emit energy, contribute to the 
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DC voltage AC voltage ) 




FIG. 2: (a)Schematic of the scattering induced by Tpi^)- 
The bandstructure of the two-channel lead L contains two 
channels p = 1,2. The time-dependent potential producing 
the scattering Tpi^) is required to be tuned so that it pro- 
duces strong transitions from channels p = 2 to p = 1 along 
the solid arrows and only weak scattering from p = 1 to p = 2 
along the dashed arrows. These transitions are separated by 
different momenta, (b) Nearly periodic grating structure to 
create the momentum selectivity needed for 7^(57). The DC 
voltages on the left end of the wire are smaller so as to con- 
fine a two-channel wire in between the solid lines, while at the 
right end they are larger so as to confine a narrower single- 
channel wire. The grating structure in the middle consists 
of a set of electrodes which are separated from each other so 
that an alternating AC voltage can be superimposed on them 
in addition to the DC voltage. The alternating AC voltage 
distorts the potential in the channel in a sinusoidal fashion as 
shown leading to an effective sinusoidal scattering potential 
between the channels p = 1 and p = 2. 



conductance in Eq. 1181 In order to achieve the required 
71^(51), which produces a conductance cr(ei?) > 1 and a 
finite pumping current, we will need to ensure that the 
scattering amplitude for the solid arrows in Fig. ^a) are 
larger than the one shown in the dashed arrows. 

One way to realize the required constraints on the scat- 
tering amplitude is to use the combination of confine- 
ment and grating potentials shown in Fig. [DJb). The 
electrodes confine the 2DEG in the middle into a quan- 
tum wire, while the grating structure allows one to apply 
a time-dependent potential, which is periodic in space 
over a certain length-scale. One can expect such a grat- 
ing structure to allow us to select specific momentum 
transfers so that the transitions along the solid arrows 



are preferred over the dashed ones. Since we will choose 
the selected momentum transfer of the grating to be of 
order Skp ~ fc_F,i — fcF,2, which is different from the back- 
scattering wave-vectors, the back-scattering rate can be 
shown to be suppressed. For simplicity, we will assume 
that the grating potential is applied with opposite po- 
larities on the two sides of the wire, while the wire is 
refiection symmetric so that the mode-wavefunctions for 
the channels p = 1 and p — 2 have opposite parities. 
Under these conditions, the applied time-dependent po- 
tential V{x,y) leading to the matrix-elements V^',p(x) 
leads to off-diagonal matrix element between the chan- 
nels p — 1 and p — 2 so that Vp>^p{x) = (1 — Sp'^p)f{x). 
Here /(x) = ^ dqdkdy<t)l ,^j^^fJ(y)V{x,y)e">''(t)2,k-q/2{y) 
and 4>p=i,2-k{y) are the possibly momentum dependent 
transverse wave-functions of the modes in the y direc- 
tion. 

The periodic nature of the grating structure in 
Fig-IHb) that generates the potential f{x), allows us to 
write f{x) = X^jl^o^ 9{^~ jc) where 5(2;) is the potential 
at the wire generated by one neighboring pair of elec- 
trodes in the grating in Fig. [UJb) with lattice constant a. 
The relevant integral over x of f{x) is written as 



dxf{x)e""' « 5(g)e*«'^(^"i)/' 



2sin(A^ga/2) 
sin {qa/2) 



(25) 

where g{q) is the fourier transform of g{x) and we have 
assumed that g{x) is such that f{x) « outside the 
interval [0,L]. Because the potential g{x) is constructed 
using electrodes with opposite signs, g{q = 0) vanishes 
at (7 = and the largest resonance (of order Ng^q = 
^)) in f{q) occurs at qa ~ 27r. The other resonances 
occuring at higher harmonics of this fundamental wave- 
length are expected to be suppressed for electrodes placed 
at a distance further than a from the wire. The width of 
the peak in f{q) is expected to be of order Sq ^ 

Performing the integrals in Eq. [24l we obtain the Flo- 
quet scattering matrices relevant to the calculation of 
o'L(e_F) in Eq. 1181 which are then written as 



. -1/2 -1/2 J./, , N 

. -1/2 -1/2 



7*(-(fco,2-A;-i,i)) 



Tf,1,1;0,2 = 
'7f,-1,1;0,2 
Tf,1,2;0,1 = 

Tf,-i,2;0,i -^^'F,V^'^^F,2'^/*(^(^o,l - ^-1,2)), (26) 



. -1/2 -1/2 J./, , N 

-IV p { VpJ /(/coa - fci,2) 

-1/2 -1/2 



where intra-band scattering processes have been assumed 
to be suppressed because of symmetry. The upper 
pair of scattering elements correspond to the solid ar- 
rows in Fig. [Ha) and the lower pair of scattering el- 
ements corresponds to the dashed arrows. It follows 
from Eq. [IS] that the conductance aLi^p) > 1, whenever 

E„' \TF,n'lM^:E = 0)\^-\TF,n'2M^,E = 0)\^ > i.e. 

the prcesses associated with the solid arrows in Fig. ^a.) 
dominate over the dashed arrows. 

The requirement that the scattering amplitudes shown 
by the solid arrows in Fig. [5] are stronger in magnitude 
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than the dashed arrows can be satisfied using the mo- 
mentum dependence of the scattering amplitude from the 
potentials Vp',p{x). This translates into the fourier trans- 
form f(q) of the x-dependence of the potential having a 
resonance of width Sq at a momentum corresponding to 
one of the solid arrows, say the upper solid arrow. To 
obtain large conductance cr(eF) > 1, the width of the 
resonance in f{q) has to ensure that the scattering asso- 
ciated with the dashed arrows in Fig. [2] are suppressed. 
The process associated with the other solid arrow need 
not be suppressed. The Fourier transform f(q) for the 
potential in Eq. [^5] has a peak of height order ^ N and 
width of order 6q ~ near qa = 2tt. By tuning a gate 
voltage, and correspondingly the mean- Fermi energy ep 
to ep — ef.Oi the Fermi wave-vectors kp^i and kp^2 of the 
incident right moving electrons in L can be tuned so that 
the condition 



kpo =^ — 



(27) 



is satisfied and all the scattering processes represented 
by the arrows in Fig. [21 in the limit of small frequency 
17 — >■ 0, can be driven by the scattering potential Vpi^p{x), 
whose Fourier transform has a peak at the wave-vector 
q = Changing the Fermi-energy ep by Sep changes 
kp.p by different rates Skpp ~ 5ep/vp,p for the chan- 
nels p = 1,2. Therefore one expects to be able to tune 
the Fermi-energy to be able to satisfy the condition in 
Eq.[27l At finite frequency fi, which is still much smaller 
than ep, the outgoing wave- vectors kn^p associated with 
the different Floquet indices n are split from kp = kp,p 
by amounts proportional to the frequency and are given 
by the relations kp^p = ko,p = + ^^^^ = - ^ 

for the channels p = 1,2. The frequency splitting of 
the wave- vector is crucial to be able achieve the regime 
where the process associate with the upper solid arrow in 
Fig-IUJa) can be made large, while keeping the dashed ar- 
rows small. The largest possible conductance, cr, should 
be obtained when the Fermi-energy ep = ep^ + Sep is 
tuned to satisfy the condition 



(fci 1 - ko 2)a = ( kp i - kp2 + ) a = 27r, (28) 



which ensures that the upper solid arrow in Fig. 
dominate the conductance in the small dq limit. In fact, 
we expect the transfer associated with the upper solid 
arrow in Fig. [SJa) to be large as long as 



n 

VPA 



27r 
a 



fSep 



n 



Vp.l 



Sep 

Vp,2 



< 



2tt 
Na 



(29) 

In addition, we require the amplitude for the processes 
associated with the dashed arrows in Fig.[21[a) to be small 



so that 

(fco.i - ki^2) = ( fcF,i - fcF,2 + — ) >— + Tr (30) 
\ vp^2 / a -'Va 

, ^ ^ , n \ 2n 27: 

(fco,! - fc_i.2) = kpA - KF,2 < 

\ Vp2j a Na 

(31) 

These conditions together with the resonance condition 
Eq. are satisfied if the frequency VI satisfies the con- 
dition 



vp,2 vpaJ Na 



(32) 



This inequality can also be viewed as a constraint on the 
total length of the grating structure 



Na> 



47ruF,ii'F,2 

ri(wF,l - Vp,2) 



(33) 



When the parameters of the system satisfy Eq. [33] and 
the Fermi-energy is tuned to satisfy Eq. [29l we expect 
the upper solid arrow in Fig. [5] to be dominant and us- 
ing Eq. [IH] it follows that (TL{ep) > 1. Using Eq. [T^ 
and Eq. 1131 leads to a conductivity (y{ep) > 1, which is 
above the Landauer limit and a finite pumping current 
ap{ep) > 0. 

The pumping current Ip{ep) at a Fermi-energy ep is 
obtained from crp(e) using Eq. 1151 In principle, eval- 
uating the pumping current Ip{ep) requires calculation 
of ap(ep) for the entire range of Fermi-energies. How- 
ever, from physical considerations, one expects a sig- 
nificant effect of the time-dependent scattering process, 
which is what leads to the pumping current to occur 
only when the Fermi-energy is tuned so that the time- 
dependent scattering is off-resonant. Therefore Ip{ep) 
should vanish when the Fermi-energy is below the range 
implied by Eq. 1291 Thus, combinging the resonance con- 
ditions and constraints in Eas. l29l3Tl Ip{ep) may be esti- 
mated by integrating (Jp{e) over a range of Fermi-energies 
ep = ep^Q + Sep around ep^Q satisfying 

Vt 2iT ^ vpi—vp2 ^ 2tt , , 
i:^<S^F— —< + 34 

wf,i Na vp^ivp^2 fF.i Na 

The corresponding range of Sep satisfies 

\Sep,nax ~ Sepmi7i\ < -77 ^'"^^''^'^ < ^- (35) 

N a vp,i — uf,2 

Given that ap^ep) < 2 in our geometry, we expect the 
pumping current Ip{e) to be bounded by and of order 
Ip ~ ^Gq for an optimally chosen set of parameters, 
where Go = 2e'^/h is the quantum of conductance. Since 
the range of Fermi-energies we have been discussing are 
relevant only for the Fermi-energy of the left-lead ep^p, 
the bound Eq. |3S] also applies to the maximum bias volt- 
age over which the linear response conductance cr{ep) is 
reasonably accurate. Thus, the bound on the pumping 



current also serves as an estimate of the excess current 
signal that one expects to see over the Landauer limit. Of 
course, the discussion of the magnitudes of the currents 
in this paragraph are rather qualitative and provide only 
order of magnitude estimates. The conductance (7{£f) 
and the pumping current in a bias window Ip can be 
calculated with quantitative accuracy by numerically in- 
tegrating the equations Eg. 1211 However, the magnitudes 
of the currents discussed in this paragraph will be quali- 
tatively verified by the numerical results presented later 
in the paper. 

Let us now discuss a choice for the various param- 
eters in our geometry that may be reasonable to re- 
alize in GaAS 2DEG experiments. We choose a lat- 
tice parameter for the grating to be a ~ 50nm. Us- 
ing Eq. [251 the choice for the lattice parameter a sets 
kF,i — kF_2 ^ = 12.56x 10^^ nm^^. Taking a nanowire 
of width ~ 35 nm^ leads to to a linear electron-density 
of TiiD ~ ^ = 0.1/ nm, which corresponds to a two- 
dimensional electron density of3.6xl0^^/ cm^ji^ Taking 
the effective mass of GaAs to be m = 0.067TOe, we obtain 
the Fermi-velocities to be vp.i = 4:Avf,2 = 2.8 x 10^ m/s. 
Choosing a frequency of ^D, ~ 40 GHz, the required 
length of the grating would be i ~ 4 fim. To achieve this 
length of grating requires N = L/a ~ 80 repetitions of 
the unit potential. 

The physics discussed in the preceeding paragraphs, 
which was based on the Born approximation, can be 
verified by numerically solving the transmission problem 
corresponding to the Hamiltonian in Eq. [21] using the 
transmission matrix technique. The results for ^{tp) — 1 
and <7p{ep) are plotted in Fig. [3| and show peaks as a 
function of energy E. These peaks correspond to the 
resonance in Eq. [211 Similar peaks occur for the other 
resonance condition corresponding to the lower solid ar- 
row in Fig. HJa). One finds that for a time-dependent 
electric field of £ ^ 1.6 x 10^ V/cm (corresponding 
to an amplitude A = max2,|/(a;)| ~ 100 /^eV), a two- 
terminal DC conductance beyond the Landauer bound 
of cr{ep) ~ 1.5 Go and a differential pumping current of 
cp(e_F) ~ 0.5 Go over a range of voltage AE ^ il. There- 
fore the total pumping current are of order / ^ 0.5ft Gq 
where Go — 2e^/h is the conductance quantum. The 
pumping current and transport current are found to be 
of order / '-^ 1.2 nA for a pumping frequency (f2/27r) ~ 40 
GHz corresponding to a voltage of 150 /^eV. 



VI. ADIABATIC PUMPING LIMIT 

In the previous sections we have shown that in spe- 
cific regimes of frequency fl, length of grating struc- 
ture L = Na and electric-field amplitude £, one can 
obtain a finite pumping current. However, we also ar- 
gued that, the pumping current must vanish in the ab- 
sence of a time-dependent perturbation. The pumping 
current Ip{ep) can be shown to vanish within the Born 
approximation, in the limit of vanishing time-dependent 




FIG. 3: The transport conductance (j{5eF + eF,o) — 1 in 
excess of the Landauer value of unity and also the pump- 
ing current Ip{SeF + eF,o) as a function of the deviation of 
the Fermi-energy, Sep = ep — eF,o, from a reference Fermi- 
energy epfi defined in Eq. [23 in the text. The deviation of 
the Fermi-energy Sep is in units of the driving frequency Q, 
while the conductance o-{ep) — 1 is in units of the conduc- 
tance quantum Go = 2e^ /h. The pumping current is in units 
of QGq. The grating structure is taken to contain A' = 80 
cycles of a sinusoidal potential and the total length is taken 
to be L = Na ~ 4/im for a frequency (n/27r) ~ 40 GHz cor- 
responding to a voltage of 150 /leV. The applied electric field, 
which is chosen to be 1.5 x 10^ V/cm (for a channel width of 
35 nm) so as to obtain the maximum conductivity. 



drive-amplitude £, by combining Eqns. [T31 HH and [HI 
An alternative approach to obtaining a static system is 
to reduce the driving frequency fl to zero, while keep- 
ing the geometric parameters of the system such as the 
length of grating L and also the electric-field amplitude 
£ fixed. Strictly speaking, this limit, which is referred to 
as the adiabatic limit, may be different from the static 
limit. However, formally the result from the previous 
section shows that Ip oc Q so that the pumping current, 
Ip, vanishes in the adibatic limit (i.e. — > 0) consistent 
with our expectations. This result must be considered a 
formal coincidence, since for fixed grating length L, the 
limit f2 — would violate the constraint in Eq. [331 so 
that one cannot simply take the limit 17 ^ in the final 
result from the previous section. 

To calculate the adiabatic limit for the current in the 
Born approximation, we expand the pumping current Ip 
in powers of fl, for small electric fields £, while the geo- 
metric parameters are held fixed. The Floquet scattering 
equations Eq. [211 are analytic in both ft and £ for fixed L 
so that the limit 17 commutes with the limit f — >■ 0. 
Therefore the leading order in £ contribution to the adi- 
abatic pumping current is calculated by substituting the 
Born scattering amplitudes from Eq. [251 into the conduc- 
tance equations Eq. [TSl and expanding to lowest order in 
the frequency fl. Using this procedure, we find that the 
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derivative of the pumping current is given by 



2 



(36) 



which vanishes as f2 — > 0. Here Skp = ^f,! — fc_F,2 and the 
fourier transform of the drive potential f{k) oc £. The 
frequency-dependence for the pumping current, Ip{€p), 
which can be calculated using Eq. I15[ is found to be 
similar to Eq. 1361 

The time-dependent driving potential in Eq. [2] can be 
written as 

{Vp^,{x) + Vlp{x)} cos m+i{VpAx)-Vl^{x)) sin int). 

(37) 

The above driving potential can be re-written as 



2[{/i?(2;)} cos int) + fi{x) sin (nt)] 



where fpix) = -^^^^Y and fi{x) - ^. 
real and imaginary parts of the function f{x). If a time- 
dependent potential with the same phase is applied to 
all the electrodes in Fig. [21 fnix) and fi{x) are propor- 
tional to each other so that effectively the spatially vary- 
ing part of the potential f{x) is real and f{k) = f*{-'k). 
This case is referred to as single-parameter pumping and 
is characterized by the constraint ^[|/(^)p — = 
0. Using Eq. 1361 it follows that the pumping current in 
the adiabatic limit for our set-up scales as Ip (x fl^. This 
scaling is qualitatively distinct from the numerical results 
presented in Fig. [31 where we found that Ip oc Q. Our 
adiabatic analysis suggests that the linear scaling of the 
pumping current /p cx 51 must cross over to an Ip cx il^ 
for sufficiently small fi, i.e. close to the adiabatic limit. 
Since the result /p cx is expected to hold only for fre- 
quencies fl ^ vp/L so that the constraint on the grating 
length L given in Eq. [531 is satisfied, one expects this 
cross-over to occur at ~ u_f/L. 

The vanishing of the linear-order in 51 pumping cur- 
rent, in the lowest-order Born approximation, is consis- 
tent with the low-frequency (i.e. — > 0) scaling of the 
pumping current from previous studies.^ For general sys- 
tems, the zero-bias current Ip in the limit of small driving 
frequency Q, can be proportional to /p cx 51 so that the 
drive pumps a fixed amount of charge per cycle of the po- 
tential. This is the reason for referring to the zero-bias 
current as a pumping current. The charge pumped per 
cycle — 



/(^)-r(^) 



(38) 



are the 



Qp = hm 

o-^o 51 



(39) 



is related to a Berry-phase of the ground-state acquired 
over a cycle of the adiabatic pump. The existence of 
a non-zero value for such a Berry phase requires the 
presence of at least two independent parameters in the 



driving potential so that the adiabatic drive can gener- 
ate a loop in parameter space. The driving potential 
discussed so far in our set-up contains only a single- 
' parameter corresponding to linearly polarized radiation. 
Therefore, based on previous results in the literature, one 
expects that the pumping charge Q p vanishes in the low- 
frequency (i.e. small 51) limit. Given the analytic nature 
of the scattering equations Eq. 1221 the pumping current 
in the case of single-parameter pumping is expected to 
scale as Ip cx 51^ as verified to o{£^) by our Born approx- 
imation result. 

The results obtained on adiabatic pumping in previ- 
ous work,— together with Eq. [36l suggest that adiabatic 
pumping with a non-vanishing linear order in 51 term 
Ip ~ Qpfl may still be possible even with significantly 
shorter structures provided a two-parameter driving po- 
tential is used. Such a potential may be implemented 
in the quantum wire set-up shown in Fig. [1] by making 
the potential f{x) complex. Physically, a complex value 
of f{x) is obtained by applying an x-dependent time- 
delay in the radiation-induced potential to the various 
electrodes used to generate the potential f{x). It follows 
from Eq. l36l that a pumping current proportional to 51 can 
be obtained by replacing the long-grating structure by a 
pair of bump potentials f{x) = e~'^'^ g{x — a) + e^'^ g{x + a) 
so that |/(fc)| = \g{k) cos {ka + ^ |/*(-A:)|. However, 
since fTp(eF) cr(ep) — 1 cx 51, the resulting deviation of 
the conductance from the Landauer limit is expected to 
be small. Therefore, the advantage in terms of conduc- 
tance excess over the Landauer limit is not expected to 
be significantly enhanced by choosing a two-parameter 
driving as opposed to a single-parameter drive. 



VII. CONCLUSION 

In this paper we have shown that the additional de- 
grees of freedom, i.e. the Floquet-index associated with 
a periodically driven system can lead to a conductance, 
which can exceed the number of channels in a quantum 
wire. For time-independent systems general arguments 
based on unitarity can be used to bound the conduc- 
tance by the number of channels in the wire. The pro- 
posed set-up consists of a two-channel wire connected to a 
single-channel wire through a time-dependent grating po- 
tential. For a gating structure of length L — Na ^ 4 fim 
in a GaAs 2DEG (seen in Fig. [3]) a peak in the conduc- 
tance beyond the Landauer limit and an oscillatory pump 
current as a function of gate voltage fj,av are calculated. 
Therefore, the Floquet-index manifests itself as an ad- 
ditional channel, which can carry a measurable current 
of the order of 1.2 nA. Furthermore, even at zero-bias 
the system carries a pumping current proportional to the 
applied frequency 51 even for a single-parameter drive. 
However, the lower bound on 51 set by Eq. [33l ensures 
that this result does not contradict the previous work 
interpreting the pumping current as the Berry-phase in 
a two-parameter space, which predicts that the pump- 
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ing current vanishes for single-parameter pumping in the 
adiabatic limit. ^ 

We thank Assa Auerbach, Dimitri Feldman, Yuval 
Gefen and Michael Moskalets for valuable discussions and 
comments. J.S. thanks the Harvard Quantum Optics 
center for support. This work was supported in part 



by NSF grant DMR 0906475. We also acknowledge the 
support from Army Research Office with funding from 
the DARPA OLE program, Harvard-MIT CUA, NSF 
Grant No. DMR-07-05472, AFOSR Quantum Simula- 
tion MURI, the ARO-MURI on Atomtronics. 



^ F. Grossmann, T. Dittrich, P. Jung, P. Hanggi, Phys. Rev. 
Lett. 67 516 (1991); M. Holthaus, Phys. Rev. Lett. 69 351 
(1992). 

^ T. H. Oosterkamp, T. Fujisawa, W. G. van der Wiel, K. 
Ishibashi, R. V. Hijman, S. Tarucha, L. P. Kouwenhoven, 
, Nature 395, 873 (1998); G. Platero, R. Aguado, Phys. 
Rep. 395, 1 (2004). 

3 D. J. Thouless, Phys. Rev. B 27, 6083 (1983) ; B. L. 
Altshuler, L. L Glazman, Science, 283, 1864 (1999) ; M. 
Switkes, C. M. Marcus, K. Campman, A. C. Gossard, Sci- 
ence 283, 1905 (1999) . 

* L. E. F. Foa Torres, Phys. Rev. B, 72, 245339 (2005). 

5 P. W. Brouwer, Phys. Rev. B 58, R 10135 (1998). 

^ M. Wagner, F. Sols, Phys. Rev. Lett. 83, 4377(1999); Y. 
Levinson, O. Entin-Wohlman, P. Wolfle, Phys. Rev. Lett 
85, 634 (2000). 

N. H. Lindner, G. Rafael and V. Galitski, Nat. Phys. 7 490 
(2011) 

® T. Kitagawa, T. Oka, A. Brataas, L. Fu, E. Demler, 

larXiv: 1 104.4636 (2011). 
^ T. Kitagawa, E. Berg, M. Rudner, E. Demler, Phys. Rev. 

B 82, 235114 (2010) 
° H. Sambe, Phys. Rev. A 7, 2203 (1973). 
^ M. A. Zudov, et. al., Phys. Rev. B 64, 201311 (2001); P. 

D. Ye, et. al., Appl. Phys. Lett. 79, 2193 (2001); R. G. 

Mani, et al.. Nature 420, 646 (2002); M. A. Zudov et al., 

Phys. Rev. Lett. 90, 046807 (2003). 



Several alternate mechanisms have been proposed to ex- 
plain these experiments. See, for example, I.M. Dimitriev, 
A. D. Mirlin, and D. G. Polyakov, Phys. Rev. B 75, 245320 
(2007); A. Auerbach and V. V. Pai, Phys. Rev. B 76, 
205318, (2007); S. A. Mikhailov, Phys. Rev. B 83, 155303 
(2011) ; and references therein. 
R. Landauer, Phil. Mag. 21, 863 (1970). 
E. N. Economou and C. M. Soukoulis, Phys. Rev. Lett. 46, 
618 (1981);i6jd. 47 972 (1981). 

D. S. Fisher, P. A. Lee, Phys. Rev. B, 23, 6851 (1981). 

A. P. Jauho, N. S. Wingreen and Y. Meir, Phys. Rev. B 50, 
5528 (1994); S. Kohler, J. Lehmann, P. Hanggi, Phys. Rep. 
406, 379 (2005); M. Moskalets, M. Buttiker, Phys. Rev. B 
66 205320 (2002); L. Arrachea, M. Moskalets, Phys. Rev. 
B 74 245322 (2006); Hernan L. Calvo, et al, Appl. Phys. 
Lett. 98, 232103 (2011); D. Martinez, R.Molina, B. Hu, 78 
045428 (2008). 

D. E. Feldman and Y. Gefen, Phys. Rev. B 67, 115337 
(2003). 

B. Kaestner, V. Kashcheyevs, S. Amakawa, M. D. Blu- 
menthal, L. Li, T. J. B. M. Janssen, G. Hein, K. Pierz, T. 
Weimann, U. Siegner, and H. W. Schumacher, Phys. Rev. 
B 77, 153301 (2008). 

G. Barak, H. Steinberg, L. N. Pheiffer, K. W. West, L. 
Glazman, F. von Oppen, A. Yacoby, Nature Physics 6, 
489 (2010). 



